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Abstract 

The abundance of Li, Be, and B isotopes in galactic cosmic rays (GCR) between E=50-200 MeV/nucleon has 
been observed by the Cosmic Ray Isotope Spectrometer (CRIS) on NASA's ACE mission since 1997 with high 
statistical accuracy. Precise observations of Li, Be, B can be used to constrain GCR propagation models. We find 
that a diffusive reacceleration model with parameters that best match CRIS results (e.g. B/C, Li/C, etc) are also 
consistent with other GCR observations. A ~ 15-20% overproduction of Li and Be in the model predictions is 
attributed to uncertainties in the production cross-section data. The latter becomes a significant limitation to the 
study of rare GCR species that are generated predominantly via spallation. 
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1. Introduction 

The enormous excess of the abundances of Li, 
Be, and B (LiBeB) relative to C, N, and O in GCRs 
by a factor of '^lO^ compared with solar system 
abundances is attributed to the fragmentation of 
primary GCR nuclei, mai nly C. N, and O, o n in- 
terstellar H and He atoms ( Freier et al. . 1959[ ). Ob- 
servations of secondary GCRs such as LiBeB thus 



provide strong constraints on propagation mod- 
els of GCRs within the Galaxy, since the produc- 
tion of these light isotopes depends on the amount 
of matter traversed during propagation. B/C ra- 
tio is often used to tune propagation parameters, 
such that the model predictions agree automati- 
cally with B/C. Li and Be are particularly interest- 
ing since their production depends not only on the 
interaction of CNO, but also on "tertiary" interac- 
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tions in the ISM (e.g. B, Be — > Li), and therefore 
may provide further restrictions on propagation 
models. In particular, one might expect to observe 
a stronger energy dependence of Li/C and Be/C 
ratios compared to the B/C ratio. In addition, es- 
tablishing the details of GCR propagation in the 
Galaxy will be beneficial for other studies, such as 
searches for signatures of exotic physics in GCR, 
spectra and origin of galactic and extragalactic 7- 
ray backgro und, nucleosynthesis, an d solar modu- 



lation (e.g.. lMoskalenko et al.l . l2005l ) 



The Cosmic Ray Isotope Sp ectrometer (CRIS) 



on the NASA/ ACE spacecraft (|Stone et al.i . ll998l) 



measures the isotopic composition of elements 
with 2<Z<30 in the energy range from ~30- 
500 MeV/nucleon since 1997 with unprecedented 
precision. High statistics allows LiBeB to be 
studied over an extended energy range (~50-200 
MeV/nucleon) for the first time. We present abso- 
lute intensities and relative elemental and isotopic 
abundances of GCR LiBeB observed by CRIS dur- 
ing near solar minimum conditions, and we com- 
pare these results with observations from previous 
instruments. We also discuss the implications of 
these observations in the context of a GCR trans- 
port model, GALPROP. GALPROP has recently 
been improved to predict LiBeB GCR isotopic 
abundances, including re-examined fragmenta- 
tion cross sections used as input to the model. 
The precisio n of the data set fro m CRIS is high 
enough (e.g.. lde Nolfoet al.l . l200ll) to require a re- 



evaluation of the contributing uncertainties in the 
model calculations, particularly the uncertainties 
in the isotopic cross-sections. 



2. Intensities and Elemental Abundances 

The absolute intensities of Li, Be, B, and C as 
measured by CRIS are shown in Fig. la (solid cir- 
cles). The data set considered for this study cov- 
ers near solar minimum conditions between 1998 
Jan. 1 and 1999 Jan. 23. Excluding periods of in- 
tense solar activity, the total tim e period corre- 



spon ds to 311 days (for details see Ide Nolfo et al 



20031 ). The uncertainties shown in Fig. la for CRIS 
observations include both statistical and system- 



atic uncertainties. The statistical uncertainties are 
typically small, ranging from 0.5% for C to 5% 
for Li. The systematic uncertainties result from 
the various correction factors within the calcula- 
tion for the absolute intensity such as uncertain- 
ties in the determination of the geometry factor 
(2%), spallation loss within the instrument (1-5%), 
and the tracking efficiency of t he Scintillating Op- 
tical Fiber Hodoscope (SOFT) (jStone et al.l . ll998l ) 
(<1% for Z>5 and up to 6% for Li). 

The elemental and isotopic ratios observed by 
CRIS are shown in Fig. lb and Fig. 2 for solar 
minimum conditions. The ratios of the spallogenic 
nuclei (B, Be, and Li) to mostly primary nuclei 
such as C are particularly important in constrain- 
ing propagation models since these ratios are sensi- 
tive to the amount of material traversed by GCRs 
from the source to detection at Earth. In addition, 
abundance ratios tend to be less sensitive to instru- 
mental uncertainties than absolute intensities. In 
Fig. lb, the experimental error bars include both 
systematic and statistical uncertainties. 

In Fig. 1, CRIS results are compared with pre- 
vious observations. A direct comparison of CRIS 
results with previous observations is complicated 
at these energies by the effects of solar modulation. 
The effect of solar modulation is approximatedby 
the spherically symmetric model of iFiskI (jl97l[ ). 
ch aracterized by the so lar modulation parameter 
(h ( Gleeson and Axfordl . [19681 . Levels of modula- 
tion are determined for a given source spectrum 
by matching post-propagation interstellar spectra 
as predicted by GALPROP to CRIS at low ener- 
gies and to IIEAO-3 at high energies. The average 
amount of solar modulation experienced by parti- 
cles observed by CRIS during near-solar minimum 
conditions corresponds to 0Ri55O MV. While the 
modulation parameter can be used to character- 
ize the amount of modulation experienced during a 
specific time period, the derived value of (j) also de- 
pends on the initial choice of the source spectrum 
and thus on the propagated interstellar spectrum 
adopted as input in the solar modulation calcula- 
tion. For instance, the solar modulation adopted 
for CRIS during solar minimum is closer to = 
~400 MV based on predictions from the Leaky Box 
Model ijPavis et al.U200ll[Nie"bur et all 



20031 ) 



CRIS observations are in good agreement with 
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Fig. 1. [a] Intensities of C, B, Be, and Li shown for two levels of solar modulation (</) = 550 (solid curve) and 800 MV (dashed 
curve)) corresponding to the epoch when the experimental data were taken, and [b] relative elemental abundances compared 
with previous observations. The experimental data is compared with the results of a propagation model GALPROP. See 
text for details. (Symbols refer to ORIS: •, Maehl et al. 1977: □, Webber et al. 1972: A, Englemann et al. 1990: ■, Orth et 
al. 1978: ♦, BufHngton et al. 1978: <), Simon et al. 1980: ♦, Buckley et al. 1994: O, Chapel and Webber 1981: V. Lezniak 
et al. 1978: A, MuUer et al. 1991: f, Dwyer et al. 1978: >, Webber et al. 1977: Garcia-Munoz et al. 1987: Duvernois 
et al. 1996: Mewaldt et al. 1981: Wiedenbeck and Greiner 1980: + , Webber et al. 2002: T, Hagen et al. 1977: *, 
Fisher et al. 1976: Q, Juliusson et al. 1974: +, Garcia-Munoz et al. 1977: O, Connell et al. 2001: X, Krombel et al. 1988: 
<, Lukasiak et al. 1999: A). 



the previous observations of Voyager 1, 2 and IMP 
7, 8, both of which were made during near solar 
modulation conditions similar to CRIS. The Voy- 
ager data were accumulated over 21 years with a 
weighted average solar modulation level of 0=450 
MV. This may not be a good representation of the 
actual modulation le vel experienced over 21 years 
( Connell et al. . 1998 h and thus complicates a com- 
parison with CRIS data. At higher energies, the 
modulation level determ ined for HEAO-3 data 
( Engl emann et al.'. 'l99Cj) corresponds to 0^800 



M V (jPavis et ah . ,200Ql . 



3. GCR Propagation Calculations 

To interpret the CRIS observations of LiBeB 
spectra and isotopic abundances, we use the dif- 
fusive reacceleration model (|Seo and PtuskinL 



19941) . "Reacceleration" is a distributed acceler- 
ation of particles due to scattering on random 
hydromagnctic waves moving at Alfven speed in 
the interstellar medium (Fermi 2nd-order mecha- 
nism). The model reproduces the peak near 1000 
MeV/nucleon in the ratios of secondary to pri- 
mary nuclei in a physically motivated way (where 
the parameters are: normalization and index of 
the diffusion coefficient, the Alfven speed, and 
the size of the galactic halo), and is consistent 
with the K- capture pare nt /daughter n uclei ratio 
2 001a: Niebur et all 



( Jones et al 



20031 ) 



While the diffusive reacceleration model suc- 
cessfully reproduces data on GCR nuclear species 



(|Moskalenko et al.l . 120021 ) . it underproduces GCR 
antiprotons by a factor of ^2 at 2 GeV, which may 
be a signature of new effects. In particular, the 
propagation of low-energy particles may be aligned 
to the magnetic field lines instead of isotropic 
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Fig. 2. Isotopic ratios of ^U/'^Li, "'B/^Be, and '^Li/^Be 
compared with previous data and with the predictions of 
GALPROP. (Data from ORIS: •, Garcia-Munoz et al. 
1977: O, Lukasiak et al. 1999: A) 



diffusion ( Moskalenko et al 1 Eooi), our local en- 
vironment (the Local Bubble) may produce a 
fresh "unpr ocessed" nuc l ei component in GCR at 
low energy ( Davis et al 1 12OOOI : iMoskalenko et al 



2003bf ). or a more intense nucleon spectra in dis- 
tant regions could yield more a ntiprotons and 
diffuse 7-rays (jStrong et all 120041 ). New accurate 
data on LiBeB spectra and isotopic abundances 
by CRIS may be used to test further th e reaccel- 



eration model ( Moskalenko et al.l . l2003a[ ). 



In our calculations we use the propagation 
model GALPROP as described in detail elsewher e 
( Strong and Moskalenkol IMoskalenko et al 



2OO2I: IStrong et all . l2004f) . The model is designed 
to perform GCR propagation calculations for nu- 
clei {Z<28), antiprotons, electrons and positrons, 
and computes 7-rays and synchrotron emission 
in the same framework. GALPROP solves the 
transport equation with a given source distribu- 
tion and boundary conditions (free escape) for all 
GCR species. This includes a Galactic wind (con- 
vection), diffusive reacceleration, energy losses, 
nuclear fragmentation, radioactive decay, and 
production of secondary particles and isotopes. 

The code includes cross-section measure- 
ments and energy dep e ndent fitting functions 
( Strong and Moskalenkol 2001 ). The nuclear re- 



action network is built using the Nuclear Data 
Sheets. The isotopic cross section database is 
built using the extensive T16 Lo s Alamos com- 



pilati on of the cross sections (jMashnik et al 
1998f l and modern nuclear cod es CEM2k and 
LAQGSM (|Mashnik et all . I2004 . The most im- 
portant isotopic production cross sections (^H, 
3H, 3He, Li, Be, B, Al, CI, Sc, Ti, V, Mn) are 
calculated usin g our fits to ma j or production 
channels (e.g.. IMoskalenko et all I2OOIL l2003bL 
Moskalenko and MashnikI 2003[ ). Other cross 
sections are calcula ted using phe n omen ological 
approximations by IWebber et al. (|l990h (code 
WNEWTR . FDR version s of 1993 and 200 3) and/or Sil- 
berberg and Tsao ( Silberberg et al.l . [l998i) (code 
YIELDX.OllOOO.FOR version of 2000) renormal- 
ized to the data where it exists. For pA inelastic 
cro ss section we adapted th e parametrization 
by Barashenkov and Polanski (code CROSEC , 



Barashenkovl 1993 : Barashenkov and PolanskH 
1994h . 

The propagation equation is solved numerically 
starting at the heaviest nucleus (i.e., ^"^Ni), com- 
puting all the resulting secondary source functions, 
and then proceeds to the nuclei with A-l. The pro- 
cedure is repeated down to A—1. To account for 
some special /3~-decay cases (e.g., ^"Be— >^°B) the 
whole loop is repeated twice. The current version 
employs a full 3-dimensional spatial grid for all 
GCR species, but for the purposes of this study the 
2D cylindrically symmetrical option is sufficient. 

For a given size of Galactic halo, matching 
the propagation model predictions to the obser- 
vations of GCR B/C both at low energy with 
CRIS data and at high e nergy with IIEAO-3 data 
l Englemann et allll990[ l determines the diffusion 
coefficient as a function of momentum and the 
Alfven speed. The halo size 4 kpc is based on 
constraints set by the r adioactive isotopes ^°Be , 
ffAL 36 ci, and ^^Mn (IStrong and Moskalenkol 
2OOII IMoskalenko et all I2OOII ). Assuming a Kol- 
mogorov spectrum of interstellar turbulence, it 
yields 6.5 x 10^^P{p/A GY)^ (cm^ s^^), where p 
is the rigidity, 6=1/3, the Alfven speed va = 35 
km s"""^, and the exact values of the reacceleration 
parameters depend mainly on the adopted cross- 
sections. The diffusion coefficient is assumed to 
be independent on the spatial coordinates. While 
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there is no a clear evidence that the diffusion co- 
efficient should be distinctly different in the disk 
and in the halo, this minimizes the number of 
parameters to be determined from the data. The 
current state of the data on radioactive isotopes 
do not warrant more a complicated approach. The 
measured isotopes to-date (^°Be, ^^Al, ^^Cl, and 
^^Mn) all have halves lives around 1 Myrs. To test 
the diffusion coefficient at different distances from 
the Sun and in the halo, one needs accurate mea- 
surements of radioactive isotopic abundances in 
CRs that include isotopes with different lifetimes 
(e.g. "'^''C). Heavy elements that have large frag- 
mentation cross sections may provide another way 
to test t he characteristics of the local interstellar 
medium (jMoskalenko et all . l2005h . 

Supernova remnants (SNR) are believed to be 
the primary sourc es of Galactic cosmic r ays. Obser- 
vations of X-rav (iKovama et al.L 
' " 120051 



1995 ) and 7-ray 



emission (jAharonian et all . I2005L l200a ) from SNR 
shocks reveal the presence of energetic particles 
thus testifying to efficient acceleration processes. 
The predicted spectrum of accelerated particles 
has power-law in rigid ity with index which may 
vary around —2.0 (e.g.. lEllison fc Cassam-Chena 



2OO5I : iBerezhko fc Volkl . |2006[ ). Such a hard injec- 
tion spectrum poses a difficulty in reconciling the 
propagated spectrum with the direct CR mea- 
surements, assuming the Kolmogorov spectrum of 
interst ellar turbulence. Th e latter i s favored by the 
data ( Woo fc Armstron"3 . 1979 : Saur fc Bieber . 
I999I ) and MHD simulations (|Verma et all . Il996l ). 
In our calculations, the injection spectrum is tuned 
to match the spectra of primary GCR nuclei for 
a given propagation model. For the reaccelera- 
tion model, the injection spectrum is taken as a 
power-law in rigidity with index 2.1 below 9 GV 
and 2.42 above 9 GV, where the injection index at 
high energies is fixed by the rigidity dependence of 
the diffusion coefficient oc p^^^ . The source abun- 
dances are tuned to match the ACE/CRIS elemen- 



tal a nd isotopic abundances (jWiedenbeck et al 



200 ll ). Generally, the injection spectrum is model 
dependent: different propagation models assume 
a different rigidity dependence of the diffusion co- 
efficient while the propagated spectra are tuned 
to the same local data. For reacceleration, this 
requires a change in the injection spectrum in- 



dex. At low energies, direct measurements of the 
interstellar spectrum are impossible due to the 
solar modulation, while indirect observations of 
the proton spectrum via pionic gamma rays by 
EGRET are not accurate enough and cannot be 
used to fix the spectrum below a few GeV. The 
spectrum of CR nucleons below ~lGGeV/nucleon 
is thus tuned to the local data assuming modula- 
tion models, which are in turn, approximate and 
also rely on the interstellar spectrum. 

With the availability of precise GCR data from 
CRIS, uncertainties in the fragmentation cross 
sections have become a significant limitation to the 
study of rare GCR species t hat are generat e d pre 



dominantly via spallat ion (jYanasak et al.l . 12001 



Moskalenko et al. ^ 200ll). A review of curr ent cross 
sections (jMoskalenko and" Mashnikl ,|200?) was un- 
dertaken for the dominant reactions involving the 
production of LiBeB as well as products decaying 
to LiBeB species (e.g., *^He, ^''^"C). For instance, 
/3^-decay of ^''Be contributes significantly to ^"^B, 
while ^Be may produce som e ^Li via electron cap- 
ture. [Sistersonet^ ( 1997 ) have surveyed partial 
^'^°Be production cross sections for fragmentation 
of O on hydrogen over a range of energies E^?>Q- 
500 MeV/nucleon. Higher energy cross-section 
measurements (i?^365-6 00 MeV/nucleon, from 
Webber et al. l ll990Lll998h for p-hCNO^LiBeB re- 



Michel et all (|l995l) 



actions have also been made, 
gives ^'^°Be production cross sections for frag- 
mentation of CNO on hydrogen over a range of 
energies i?^800-2600 McV/nucleon. In addition, 
measurements for He-t-CN-^BeB re actions at 600 
MeV/nucleon dWebber et all. Il990l) an d a + a^ 



19971 ) are 



Li, ■^Be at 60-160 MeV ([Mercer et al 
available. 

Li, Be, and B in cosmic rays are mainly produced 
from the interaction of CNO nuclei with interstel- 
lar hydrogen. Some of the dominant reactions, such 
as -> 6Li, 7^9>iOBe, "B, ^ i^q ^ 
^Li, ^'^"Be, ^'-'^^B, are constrained by the data. The 
cross sections for the production of boron are bet- 
ter known than for the lighter isotopes of beryllium 
and lithium, since the main contribution to boron 
production is from CNO and ^^B — > B. Thus the 
model parameters are rather tightly constrained 
by the measured B/C ratio at low and high ener- 
gies. However, some channels of Be and Li show a 
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spread in the data between 10-20% or greater. For 
example, ^Be production by CNO nuclei shows a 
spread by a factor of ~2 at energies below a few 
100 MeV/nucleon. Some important channels are 
represented by one or few data points in a narrow 
energy range (e.g., "B, C, O ^ ^^B, "B ^Be), 
while some are not measured at all (e.g., ^"^B, N 
^Be, B, Be ^Be). The production of Li iso- 
topes is particularly poorly measured. There are 
only a few data points for Li production by CNO 
nuclei. The data on Li isotopic production by ^Li, 
Be, B nuclei are absent, while a contribution from 
spallation of these isotope s may exceed ~35%. See 
Moskalenko and MashnikI (j2003r ) for more details. 

When there is no data, the production cross 
sections are calculated using the semi-empirical 
parametrizations. This translates to about a 15% 
overall uncertainty in the isotopic production 
(~20% for Li), given that typical cross section un- 
certainties of semi-empirical parametrizations are 
of the order of -50%. 

The solid curves in Figs. 1 and 2 correspond 
to GALPROP predictions modulated assuming a 
modulation parameter of 550 MV, corresponding 
to the time period covered in this study (Jan. 1, 
1998 to Jan. 23, 1999). The dashed curves in Fig. la 
refer to modulated GALPROP predictions with 
(p—SOO MV, appropriate for the modulation level 
during HEAO-3 observations. The model is tuned 
to match the observations of C and B intensities 
and the B/C ratio from ORIS (Fig. la/b). How- 
ever, Be and Li intensities and the ratios Li/C and 
Be/C are ~10"20% lower than predicted by the 
model, although still consistent within the uncer- 
tainty determined for the propagation calculation 
(shown as the hatched region) . The isotopic ratios 
shown in Fig. 2 are also in good agreement with 
the model predictions. The model uncertainties are 
mostly due to uncertainties in the cross sections of 
tertiary interactions such as the production of Be 
and Li from isotopes of Li, Be, and B as discussed 
above. Further work is required to pin down these 
uncertainties. 



4. Conclusion 

We have made GCR LiBeB abundance mea- 
surements using the CRIS instrument during 
near solar minimum conditions in 1998-1999. The 
isotopic ratios of LiBeB are in agreement with 
previous measurements of GCR light isotopes, 
particularly for those experiments acquiring data 
during periods of solar modulation levels similar 
to the modulation levels in 1998-1999 covered by 
CRIS. The reacceleration model gives a satisfac- 
tory prediction for GCR primary and secondary 
species with Z>3, and shows good agreement with 
relative isotopic abunda nces ( Jones et al.l . l200l"bt 
iMoskalenko et all . 12002 ). In some cases, notably 
(p, a)-l-(Li,Be,B)^Be,Li, a lack of cross-section 
measurements limits the reliability of model pre- 
dictions. Understanding the subtle differences 
between the model predictions and experimental 
data now hinges on more precise cross section 
measurements, especially for reactions involving 
the production of the light isotopes of LiBeB. 
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